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The motion of polymers with inhomogeneous structure through nanopores is dis- 
cussed theoretically. Specifically, we consider the translocation dynamics of polymers 
consisting of double-stranded and single-stranded blocks. Since only the single- 
stranded chain can go through the nanopore the double-stranded segment has to 
unzip before the translocation. Utilizing a simple analytical model, translocation 
times are calculated explicitly for different polymer orientations, i.e., when the single- 
stranded block enters the pore first and when the double-stranded segment is a lead- 
ing one. The dependence of the translocation dynamics on external fields, energy 
of interaction in the double-stranded segment, size of the polymer and the fraction 
of double-stranded monomers is analyzed. It is found that the order of entrance 
into the pore has a significant effect on the translocation dynamics. The theoretical 
results are discussed using free-energy landscape arguments. 
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I. INTRODUCTION 

Translocation of polymer molecules through nanopores plays a critical role in nu- 
merous chemical and biological processes.—*^ Examples include the motion of DNA and 
RNA molecules across nuclear pores, injection of viral DNA into cells, gene swapping 
and protein transport through membranes. Because of its importance for understand- 
ing fundamental processes in biology^ and potential technological applications'^ the poly- 
mer translocation has become a subject of intensive experimental,— 
theoretical^'^i^^i^«i^9i^°i^^i^^i^-'^i^^i^-^i^6 numerical^i^^i^^i^^i-'^Qi-'^^i-'^^i^''^ studies. 

In in vitro experiments the polymer chain threads through the single protein channel, 
like the bacterial a-hemolysin, that is positioned in a 

membrane .^'S'£'2»S<miLi2ii^ For the 

polymer moving through the pore the number of available configurations decreases, creating 
an effective entropic barrier for translocation. To overcome this barrier the polymers, which 
are typically DNA or RNA molecules with charged bases, are driven across the membrane 
with the help of external electric fields. When the polymer is in the pore the membrane 
current, produced mainly by small ions present in the system, decreases significantly. The 
translocation properties of every polymer molecule is specified by measured translocation 
times and blockade currents. Because the narrowest part of a-hemolysin channel is less 
than 2nm in diameter , ^i^i'^'i^i^i^^i'^^i'^^i'^^i'^^ only single-stranded DNA or RNA molecules can 
pass through the channel. In addition, the translocation for not very long DNA and RNA 
molecules is typically much faster than other biologically important processes, such as inter- 
and intra-molecular interactions.— Explicit measurementsi^i showed that for the range of 
voltages from 30 to 150 mV the translocation time for the single-stranded DNA with 26 
bases is less than 0.01 of the time needed to unzip a hairpin with 10 bases under the 
same conditions. These observations allow to utilize the nanopore experiments to test the 
structure, chemical states and interactions of single polymer molecules. For example, in Ref.— 
single-molecule dynamics of unzipping of double-stranded DNA has been investigated by 
pulling one of the strands through a-hemolysin pore at different voltages and temperatures. 
Other recent experiments^^ also show an orientation asymmetry in the translocation of 
single-stranded DNA molecules, which is due to complex interactions of DNA nucleotides 
with the protein nanopore. 

Theoretically, the translocation dynamics is generally viewed as an effective one- 
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dimensional diffusion in a free-energy landscape that is biased by en external driving 
field . The translocation times are associated with mean first-passage 
times that can be calculated from the free-energy potentials . ^^'^^ Most of the theoretical 
works consider only homogeneous polymers, although several nanopore experiments^'^^i'^^^ 
indicate that in the real biological systems inhomogeneities in the structure and interac- 
tions of translocating polymers with other molecules might have a significant effect on the 
overall dynamics. Muthukumar— theoretically discussed some effects of the sequence on 
DNA translocation through nanopores by considering the motion of diblock copolymers. It 
was found that for the weak driving forces the passage time depends on which polymer end 
begins the translocation. However no explanation for this phenomenon has been given. The 
effect of sequence heterogeneity in polymer translocation was also studied in Ref.— where it 
was shown that the heterogeneity might lead to anomalous dynamics at some conditions. 

In this paper we investigate the translocation dynamics of inhomogeneous polymers using 
simple analytical model and free-energy landscape arguments. Specifically, we consider the 
motion of polymers consisting of two blocks. One is a single-stranded chain and another 
one is a double stranded segment. Translocation times at different external conditions are 
calculated explicitly for two cases: 1) when the single-stranded part moves through the 
nanopore first and then the double-stranded block unzips and translocates; and 2) when the 
double-stranded end comes first to the pore, which means that this segment starts to unzip 
and then one of the unwounded chains enters the pore followed by the translocation of the 
single-stranded piece. 

The paper is organized as follows. The description of the model and our theoretical 
method are presented in Sec. II, while the results and predictions are discussed in Sec. 
III. Section IV summarizes and concludes our theoretical analysis. Some details of the 
mathematical calculations are given in Appendix. 

II. OUR MODEL 

Let us consider a translocation through the pore of the inhomogeneous polymer molecule 
with monomers as shown in Fig. 1. Each monomer has a charge q. One part of the 
molecule is a single-stranded segment, while the other part is made of double-stranded 
segment. The composition of the polymer is specified by a parameter a (0 < a < 1), i.e., 
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Na monomers belong to the double-stranded block with a bond energy per monomer e. We 
assume that the diameter of the nanopore is so small that only the single-stranded polymer 
chain can go through the channel. If the polymer comes to the entrance of the pore with 
the double-stranded segment first, then it starts to unzip and only the longest chain goes 
through the nanopore - see Fig. 1. It is also assumed that the length of the nanopore is 
small so that no more than one monomer can be found inside. It is a good approximation for 
the polymers whose contour length is much larger than the nanopore length. Note, however, 
that the length of the nanopore is an important factor in the translocation dynamics of 
real polymers . ^'^^'^^ In addition, we discuss only the translocation process after the polymer 
entered into the pore, and it is assumed that the molecule that entered into the pore will 
not exit back. The problem of searching for the nanopore is beyond the scope of this work. 

To overcome the entropic barrier the polymer is driven across the channel by an external 
electric field of strength V. In the model the hydrodynamic force is neglected since theoreti- 
cal calculations show that it is smaller than the electric force in the in vitro experiments.— "^i 
We define Pk{t) as a probability to find k monomers already translocated at time t. The 
dynamics of the system at all times is governed by a Master equation. 



where Uk (wk) is the rate to move the polymer molecule by one monomer in the forward 
(backward) direction. The rates are related by the detailed balance condition. 



where Fk is a free energy of the polymer with k monomers already moved across the nanopore 
and (3 = l/Zc^T. It is also assumed that if the free-energy difference of moving the monomer 
forward or backward is zero, i.e. AF^ = for all k, the forward and backward rates are the 
same. 

The free energy of the translocating polymer, shown in Fig. 1 has three terms, namely, 
entropic, electrostatic and the hybridization energy. 



Uk-lPk-l{t) + Wk+lPk+l{t) - {Uk + Wk)Pk{t) 



(1) 



dt 




(2) 



Wk+l 



Fk Fk^entr ~l~ Fk^elec ~\~ Fk^hyb- 



(3) 



However for realistic values of external fields and unzipping energies the entropic contribution 
can be neglected.— >2i Specific values of the free energy depend on the entrance order into the 



5 



nanopore. For situations when the single-stranded segment is entering first the free energy 
is given by 

-kqV-eNa, for k < N{1 - a); 
-k{qV-e)-eN, for k> N{l-a), 
where q is an effective charge per monomer, and e is a hybridization energy per monomer. If 
the double-stranded segment is first at the entrance into the nanopore, then the free energy 
of the system is different. It also a time- dependent function. Before full unzipping of the 
double-stranded segment we have 

Ffc = -k{qV -e)- sNa, for k < Na. (5) 

We assume that the process of full unzipping is irreversible, and the free-energy potential 
changes after that moment. It becomes equal to Fk = —kqV for all values of k. The analysis 
of Eqs. (jH) and (0) indicate that at e = qV the free-energy surfaces change the behavior. 
As we will see later, it has a strong effect on the translocation dynamics. 

Dynamic properties of translocating polymers depend on the forward and backward rates. 
However, the detailed balance condition Q specifies only the ratio of the rates. In this work 
for simplicity we assume that the free energy change influences only the forward rates. The 
rates are given by 

Uk = w'exp[PqV], Wk = w', (6) 

for moving the single-stranded segment across the pore, and for transporting the single chain 
after unzipping the double-stranded segment, 

Uk = wexp[(3{qV - e)], Wk = w. (7) 

More general cases when the free energy change can be distributed between the forward and 
backward rates can be easily accommodated in our theoretical approach, however, quali- 
tatively it does not change the results. In Eqs. ()(i|7j) parameters w and w' are attempts 
frequencies, and they describe the motion in the absence of external forces and hybridiza- 
tion. Note, that in our simplified model the asymmetric interactions between the pore 
and the polymer~ are not considered, although they can be easily incorporated. However, 
since it is found that the unzipping controls the translocationi^, we believe that asymmetric 
interactions are less important in our case. 
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In this theoretical approach both forward and backward transitions in the translocation 
process are taken into account. These transitions are easy to understand for the motion of the 
single-stranded segment. However, there is some ambiguity for the double-stranded region. 
We assume here that the forward step, when the doubles-stranded region moves, corresponds 
to unzipping that is followed by the forward translocation of one monomer. Similarly, the 
backward step describes the reverse translocation followed by zipping. Obviously, other 
processes can take place in the system but the model describes the most probable events of 
translocation of diblock copolymer. Recent theoretical calculations show that at low external 
fields the backward transitions out of the pore might become an important factor.^^ 

The dynamics of the polymer moving through the nanopore is described by its translo- 
cation time. We define ri as the translocation time of the polymer with the single-stranded 
segment entering first into the channel. Correspondingly, T2 is the translocation time when 
the single-stranded segment moves into the pore after unzipping and translocation of the 
double-stranded block (see Fig. 1). Translocation times are calculated exactly for given pa- 
rameters using the standard expression for mean first-passage times.-^^^ Then to compare 
two different translocation mechanisms we introduce a dimensionless ratio of passage times, 

r = n/T2. (8) 

III. RESULTS AND DISCUSSION 

The passage times are computed exactly for the different sets of parameters. Explicit 
expressions are quite bulky and presented in Appendix. The analysis of formulas ()A2|A4|) 
suggests the translocation times can be written as a sum of three terms, 

Ti = t{s) + T{d) + Ti{corr), i = 1,2, (9) 

where r(s) is the time to move only the single-stranded segment of the polymer and T{d) 
corresponds to the time to thread only the segment that was double-stranded before the 
translocation. Obviously, these two contributions are independent of the orientation of the 
translocating polymer. The last term in Eq. Q is a correlation time that reflects backward 
trajectories when some part of the polymer already unzipped in the previous steps might zip 
again. It depends on the orientation of the polymer in the nanopore, and it plays a critical 
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role in the translocation dynamics. 

Translocation dynamics depends on the external voltage V, hybridization energy e, the 
size N of the polymer molecule and the fraction a of the double-stranded monomers. Exact 
analytical formulas allow us to consider these effects in detail. In the limit of high external 
voltages the translocation times in both orientations become equal, and they exponentially 
approach zero as follows, 

n = r2 ~ ^^il^ expi-/3qV) + — exp[/?(£ - qV)] ^ 0. (10) 

Because of the very strong field, in this regime the probability of backward transitions is very 
low, and the correlation time is essentially equal to zero. This leads to the dynamics that 
does not distinguish between the orientation of the translocating polymer. The situation is 
very different for the weak external forces, where the explicit formulas ()A2|A4|1 yield, 
n = —[N^l - af - N{1 - a)] + ^{1 - a)[l - eMPNae)] 



and 



T, = ^[{l-a')N'-Nil-a)] + A, (12) 



with an auxiliary function A defined as 



Na , Nil -a) Naexp(2Pe) 

A = exp(/3e) + ^- — ^+ ^ 



w w' w[l — exp{pe)\ 

exp[pe{Na + 1)] exp{2Pe) exp[Pe{Na + 1)] . . 

~ w[l - exp{f3e)] " w[l - exp{(36)f ^ w[l - exp{f3e)]'^ ' ^ ^ 

Passage times to move a polymer molecule across the channel from two opposite configu- 
rations differ at low external fields. These results are illustrated in Fig. 2 where the ratio 
of translocation times r is plotted as a function of external voltage. At qV/ksT ^ 1 the 
correlation terms in the translocation times differ significantly for the two orientations. This 
can be understood by considering the situation when one block is already fully translocated. 
If the translocated block is the single-stranded one, then the forward motion is energetically 
unfavored since it require to unzip double-stranded monomers. As a result, there will be a 
lot of backward transitions. In the opposite orientation, when the double-stranded region 
already passed the channel, the situation is different. Here the probability of forward and 
backward steps is equal and not so many backward trajectories are observed. As a result 
we have ti > T2, and it was observed that this inequality generally holds for all parameter's 
space. 
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Varying hybridization energy has also a strong effect on the translocation dynamics. 
When e = two orientations are indistinguishable (for w = w') and, as expected, r = 
ti/t2 = 1. However, for strong interactions in the double-stranded segment it is important 
which polymer segment enters the nanopore first. From Eqs. ()A2|A4jl it can be shown that 
for large e 

exp[(3{Na + 2){e- qV)] exp[/3(e - 2qV)] exp[(3Na{e - qV)] {exp[-(3N{l - a)qV] - 1) 



and 



w[l -exp(/?(e - g\/)]2 w'[l - exp{(3{e - qV)][l - exp{-f3qV)] 

(14) 



exp[f3{Na + 2){e-qV)] 



w[l — exp[/3{e — qV)] 
Then the ratio of translocation times for e ^ 1 is given by 

w{exp{-i3qV) - exp[-/3(iV(l - a) - l)qV]) 

w'[l - exp{-PqV)] ■ ^ ^ 

The dependence of r on hybridization energy is plotted in Fig. 3. It can be seen that 
strong intra-molecular interactions significantly slow down the translocation of polymers 
with single-stranded segment entering first. This is due to the fact that the slow unzipping 
allows the single-stranded polymer chain to make many forward and backward threading 
transitions. It leads to ti being typically much larger than T2 at these conditions. This effect 
is even stronger for weaker external voltages. 

An important parameter in the motion of polymers across the channels is the molecular 
size A^ .^i2ii22 Q^j^ theoretical calculations for diblock copolymers with — oo show that 
the translocation depends on the relative values of hybridization energy and electrostatic 
energy. From Eqs. ()A2IA4|) for e > qV the passage times are equal to 

_ exp[p{Na + 2){e-qV)] _ exp[/3{e{Na + 1) - qV{Na + 2))] 

w[l - exp[p{e - qV)Y w'[l - exp[p{e - qV)][l - exp{-PqV)y ^ ' 

and 

_ exp[(3{Na + 2){e-qV)] 

w[l - exp[f3{e - qV)]^ ' ^ ^ 

It is interesting to note that at these conditions both ti and T2 depend on A^ only through 
Na. This means that the translocation times are completely determined by the size of the 
double-stranded segment and they do not depend on the size of the single-stranded block. 
The translocation dynamics changes dramatically for large-size polymers when the electric 
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field starts to dominate, i.e., when e < qV . In this case in both orientations the polymer 
moves through the nanopore with equal times, 

iV(l -a)exp(-/3gV) Nae^^[l3{£ - qV)] A^(l - a) exp(-2/?gr) 

w' w w'[l — exp{—PqV) 

Naexp[2(3ie-qV)] 

w[l - exp[f3{e - qV)y ^ ^ 

This is exactly the sum of the passage times for two polymer blocks separately, i.e., the 
correlation term in the translocation time [see Eq. (jH))] is negligible at these conditions 
when the electric field dominates the passage dynamics. Thus we obtain for ^ 1, 

1, e <qV; 

(20) 

1 I w exp[/3(£-2qy)] ieMl3{e-qV)] - 1) ^ ^ „T/ ^ ^ 

«)'cxp[2/3(e-gV)][l-cxp(-/3gV)] ' ^ ^ V " • 

The fact that e = qV specifies a dynamic phase boundary could also be understood from 
free-energy potentials - see Eqs. Q and (jSJ. The ratio of translocation times as a function 
of polymer size is shown at Fig. 4. This situation is an example of a dynamic phase 
transition (A^ oo) when changing the electrostatic energy contribution (by modifying V 
and/or the effective charge q) leads to a transition from the orient at ion- dependent transport 
to the dynamics that does not discriminate between the polymer orientations. For realistic 
systems e ~ 2-3 ksT, and for typical translocation voltages V = 100 mV and q = le 
we obtain at room temperature qV ~ 4 ksT. However, the charge on the monomers 
can be significantly lower due to screening." This suggests that the predicted dynamic 
transition might be observed in the nanopore experiments by modifying the voltage or the 
ionic strength (to change the charge). The main observable feature in this transition is a 
width of the passage-time distributions. For orientation-dependent dynamics the width is 
wide because of possible overlapping of two distinct distributions. At the same time, the 
width narrows for the dynamics that does not depend on the polymer orientation. It can 
be suggested that this result possibly can be used as an experimental method of measuring 
effective charge q at each monomer. Furthermore, this phenomenon could be important in 
biological processes. 

Another parameter that can influence the translocation dynamics is the fraction a of the 
double-stranded monomers. The situations when a = and a = 1 are trivial since they 
correspond to the case of only single-stranded or double-stranded polymers, for which there 
is no orientation dependence. The general case < a < 1 is more complex with a non- 



10 



monotonous behavior, as can be seen in Fig. 5. The maximal discrimination between two 
different polymer orientations depends on the hybridization energy and the external voltage. 

Our theoretical calculations of the translocation of diblock copolymers predict several 
unusual orientation-dependent phenomena. This dependence can be explained using free- 
energy landscape arguments. The polymer molecule before entering the nanopore has the 
same free energy in both configurations. Similarly, after moving completely through the 
channel there is no difference in free energy of the polymer in both orientations. However, 
when only the part of the molecule translocated the free energy of the system depends on 
what segment of the polymer entered first into the pore. This indicates that, although the 
initial and final points on the frcc-cncrgy landscape arc the same, polymers in different 
orientations follow different free-energy pathways, leading to orientation-dependent phenom- 
ena described above. This is similar to a chemical reaction with and without catalyst. The 
catalyzed reaction is faster because it follows a different free-energy pathway. 

IV. SUMMARY AND CONCLUSIONS 

We developed a simple theoretical model for the translocation through the nanopore 
polymer molecules consisting of single-stranded segment and double-stranded segment. Us- 
ing exact analytical results for the mean first-passage times we show that for some sets of 
parameters the translocation dynamics of inhomogeneous polymers depends on the orienta- 
tion of the molecule at the entrance. The translocation for two orientations is different al 
low external forces, while it disappears at large voltages. Strong hybridization energies also 
lead to orientation-dependent dynamics. The effect is stronger for large polymer sizes. In 
addition, we also observed an interesting dynamic phase transition at large by changing 
the relative contributions of external forces and hybridization energy. Our theoretical results 
are discussed in terms of the free-energy landscape arguments. 

The presented theoretical model of translocation is rather oversimphfied. In order to 
have more reahstic description several factors should be included in the analysis. These 
factors include entropic contributions to free energy that are important at weak external 
forces, the effect of sequence dependence in rates and hybridization energies, and the effect 
of folding and unfolding before and after translocation. It will be also very interesting to 
fully investigate the observed dynamic transition. 
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Appendix 

Translocation times can be calculated using the well-known exact formula for mean first 
passage times for a one- dimensional asymmetric discrete random walk.— For the case when 
the double-stranded block is first at the entrance into the pore it is given by 

Na-l 



r, = ^ expiPie - qV)] + i^^il^ exp(-/?g\/) + ^^P[/^(^ " g^)] y (ATa - k) exp[pk{e - qV)] 

k=l 

+ ^^^PtM)iV(l-a)Vexp(-/5A:gT/) 



w 

k=l 

(n N{l-a)-l 

^ exiHJ^gv^ ^ [iV(l - a) - A;] exp[-/5(Ara + (Al) 

^ k=l 

Direct summation then yields the following final expression, 

T2 = exp[/3(e - gV)] + ^ — ^-exp{-PqV)+ '^^ 



w w' w{l — exp[/3{e — qV)]) 

exp[l3{Na + l){e - qV)] _ exp[2f3{e - qV)] exp[/3(A^a + l)(e - qV)] 
w[l - exp[(3{e - qV)] ~ w[l - exp[(3{e - qV)]"^ ^ w[l - exp[(3{e - qV)]"^ 
exp{-2f3qV) exp[-l3{N + l)qV] exp[-(3{Na + 2)qV] 



2 



+ ^^w'[l - exp{-pqV)] w'[l - exp{-PqV)] w'[l - exp{-PqV)] 

exp[-(3{N + l)qV] 

w'[l-exp{-(3qV)Y' ^ ' 



Similar calculations for the single-stranded segment translocating first produce 



n = ^ expme - qV)] + ^^^-"^ expi-PqV) + ^"P^^^^'^^^^ V (Na - k) exp[/5fc(^ - <lV)] 

k=l 

in T/N N{l-a)-l Na N{l-a) 

^exjp^jjqv^ ^ [Ar(l-a)-fc]exp(-/3fcgV-) + -^ ^ exp[/3(te-(^ + %V)],(A3) 



k=l k=l i=l 
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that can be written in the following final form after summation of all terms, 

Na iV(l-a) , ^ exp[{Na + 2)(3{e-qV)] 

w w' w[l — exp[/3(£ — qV)\^ 

^ Na cxi)[2,](--(/V)] exi)[2,](-- 

w 1 — e-xjp[l3{e — qV)] wfl — exp[/9(£ — gl^)]^ 

+ exp -/3A^ 1 - a)qV]^ + 

w'[l — eyiY>[—pqV)\^ w' l — exp[—(jqV) 

exp{-2pqV) exp[p{e - 2qV)] 



+ 



w'[l-exp{-PqV)]^ w'[l-exp[p{e-qV)]][l-exp{-PqV)] 

exp[pNa{e - qV) - (3N{1 - a)qV] exp[/3{e - 2qV)] cxp[-(3N{l - a)qV] exp[l3{e - 2qV)] 



w'[l - exp[i3{e - qV)]][l - exp{-(3qV)] w'[l - exp[(3{e - qV)]][l - exp{-f3qV)] 

exp[f3Na{e - qV)] exp[f3{e - 2qV)] ^^^^ 



w'[l - exp[P{e - qV)]\[l - exp{-PqV)] ' 
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Figure Captions: 

Fig. 1. Schematic view of moving a diblock copolymer, made of single-stranded and 
double-stranded segments, across the nanopore. Two different polymer orientations are 
shown. Only single-stranded chain can travel through the nanopore. Na is number of 
monomers in the double-stranded region, a) Translocation with the single-stranded segment 
entering the nanopore first, b) Translocation with the double-stranded segment moving first 
and unzipping. 

Fig. 2. Ratio of translocation times as a function of external electric field. The hybridiza- 
tion is £ = 2 ksT. Sohd lines describe polymers with = 30 monomers, while dashed lines 
are for polymers with N — monomers. In each set upper curves correspond to a = 0.25, 
and lower curves are for a = 0.75. For all calculations it is assumed that w = w' . 

Fig. 3. Ratio of translocation times as a function of hybridization energy e. Polymer 
molecules with iV = 30 monomers is considered. For solid lines a = 0.75, while for dashed 
lines a — 0.25. In each set upper curves correspond to qV/kBT — 0.05, and lower curves are 
for qV/ksT = 0.1. For all calculations it is assumed that w — w'. 

Fig. 4. Ratio of translocation times as a function of polymer size N. The hybridization is 
e — 2 ksT. Solid lines describe the translocation with qV/kBT = 0.1, and for dashed lines 
qV/ksT — 0.05. In each set upper curves correspond to a = 0.25, middle curves correspond 
to a = 0.5, and lower curves are for a = 0.75. For all calculations it is assumed that w = w'. 

Fig. 5. Ratio of translocation times as a function of the relative size of the double- 
stranded segment. For calculations it is assumed that e = 2 ksT, qV/ksT = 0.1 and 
w = w' . The upper curve describes the polymer with = 60 monomers, while the lower 
curve corresponds to A?" = 30. 
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Figure 2. Kotsev and Kolomeisky 
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Figure 3. Kotsev and Kolomeisky 
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Figure 4. Kotsev and Kolomeisky 
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Figure 5. Kotsev and Kolomeisky 



